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Abstract. Five years of EROS data towards the Small Magellanic Cloud have been searched for gravitational 
microlensing events, using a new, more accurate method to assess the impact of stellar blending on the efficiency. 
Four long-duration candidates have been found which, if they are microlensing events, hint at a non-halo popula- 
tion of lenses. Combined with results from other EROS observation programs, this analysis yields strong limits on 
the amount of Galactic dark matter made of compact objects. Less than 25% of a standard halo can be composed 
of objects with a mass between 2 x 10"*^ Mq and 1 Mq at the 95% C.L. 

Key words. Galaxy; halo - Galaxy; kinematics and dynamics - Galaxy: stellar content - Magellanic Clouds - dark 
matter - gravitational lensing 
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properties of the clouds could account for differences if the 
events are due to self-lensing. 

The EROS 2 experiment has been surveying the SMC 
since 1996. We present here an analysis of five years of 
data accumulated towards this target. 



2. Experimental setup and SMC observations 

The telescope, camera, teles cope operations and data re- 
ductio n are as described in IPalanaue-Delabrouille et alJ 
l)l998|) . hereafter SM98, and references therein. Since July 
1996, 10 one-square-degree fields have been monitored to- 
wards the SMC. The first 5 years of data from 8.6 square 
degrees spread over these 10 fields have been analyzed. 
The data set contains about 5.2 million pairs of light 
curves in two wide pass bands called thereafter "red" and 
"blue", sampled on average once every 2.5 days from end- 
April to mid-March when the SMC is visible. About 400- 
500 images of each field were taken, with exposure times 
ranging from 5 minutes in the center of the SMC to 15 
minutes in the outermost region. 



3. Data analysis 

The analysis of the five-year data set is similar to those 
of the first v ear (sm98) and of the first two years 
I Afonso et al of data where details can be found. 

A first set of cuts is designed to select light curves 
that exhibit a single significant fluctuation. This already 
excludes 99.8% of the stars, most of which exhibit flat 
light curves. 

A second set of cuts rejects stars that are a priori 
likely to be variable, based upon their position in a color- 
magnitude diagram: bright blue stars in the upper main 
sequence with low-amplitude variations and stars much 
brighter and redder than those of the red clump arc re- 
jected. These stars are contained in two thinly populated 
regions of the color-magnitude diagram, so rejecting them 
does not reduce much the number of light curves on which 
we can search for microlensing. Stars exhibiting a strong 
correlation between the red and the blue light curves out- 
side the period containing the main fluctuation are also 
rejected. Supernovae are rejected by the comparison of 
the rising and falling times of the variation. Half of the 
light curves surviving the flrst set of cuts are rejected at 
this stage, leaving about 5000 light curves. 

A third set of cuts improves the signal-to-noise ratio of 
the set of selected candidates by comparing the measure- 
ments with the best-fit point-source point-lens microlens- 
ing light curve. These criteria are sufficiently loose not to 
reject light curves affected by blending, parallax and most 
cases of binary lenses or sources. The remaining sample 
consists of about 30 light curves exhibiting a clean and 
unique fluctuation with a smooth time variation. 

A fourth set of cuts constrains the time coverage of 
the event, excluding in particular all events with Einstein 
radius crossing times Ie exceeding 1200 days. 



The final criterion is as follows: if blending significantly 
improves the fit, then the fitted blending should be phys- 
ical, i.e. indeed correspond to the amplification of only a 
fraction of the total fiux recovered. 

To limit contamination of the set of selected candidates 
by spurious events, the levels of the above cuts are set such 
that most variable stars are rejected by at least two of 
them. The tuning of each cut and the estimation of the ef- 
ficiency of the analysis is done with Monte Carlo simulated 
light curves, as described in SM98. Microlensing parame- 
ters are drawn uniformly in the following intervals: time 
of maximum magnification to within the observing period 
±300 days, impact parameter normalized to the Einstein 
radius uq G [0, 1.5] and time-scale G [0, 1200] days. 

Four candidates (SMC-1 to SMC-4) pass all the cuts. 
Their light curves are shown in figure ^ their coordinates 
are given in table ^ and the microlensing fit parameters 
(point-source, point-lens and no-blending) in table |2| 
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Fig. 1. Red light curves of the microlensing candidates 
detected towards the SMC. 

Candidate SMC-1 is the one already present in the pre- 
vious analyses (one year and two years of SMC data) done 
by EROS. Further analysis of its light curve indicated that 
it was probably a self-lensing event (sm98). It was also de- 
tected as a possible microlensing candidate bv the online 
trigger system of the MACHO experiment l|Alcock et al.L 
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Table 1. Coordinates (J2000) and baseline photometry of 
the candidates. 





a 


S 


Vj 


Ic 


SMC-1 


01:00:05.73 


-72:15:02.33 


17.9 


17.9 


SMC-2 


00:48:20.16 


-74:12:35.31 


20.4 


19.5 


SMC-3 


00:49:38.68 


-74:14:06.77 


20.2 


19.3 


SMC-4 


01:05:36.56 


-72:15:02.39 


19.9 


19.4 



Table 2. Results of microlensing fits to the SMC candi- 
dates, uo is the impact parameter and the Einstein 
radius crossing time in days. 





uo 




xVd.o.f. 


SMC-1 


0.52 


101 


1026/706 


SMC-2 


0.82 


390 


1526/990 


SMC-3 


0.66 


612 


2316/968 


SMC-4 


0.80 


243 


2923/951 



[19973). Candidates SMC-2, SMC-3 and SMC-4 have very 
long durations, with few points outside the amplified re- 
gion of the light curve. The x^/d.o.f. of candidate SMC-4 
is large, yet far from the cut level (which is set at 5.0 in 
the peak region). All three candidates, however, exhibit 
features in their baseline that are reminiscent of variable 
stars. These features appear clearly on the light curves 
plotted in 25-day bins (see figure [21). To allow a direct 
comparison of the red and blue fluxes, all fluxes have been 
normalized so that they would have a mean of zero and 
a variance of unity. Candidates SMC-2 and SMC-4 clearly 
look like irregular variables, while SMC-3 has an irregular 
light curve and seems to be rising at the end. These can- 
didates will need to be monitored for several more years 
to be either confirmed but most probably ruled out as mi- 
crolensing candidates. They are not in the ~3 square de- 
grees of the SMC on which the MACHO project took data. 

The efficiency of the analysis for events with (Monte- 
Carlo) Mo < 1.5 normalized to an observing period Tobs of 
five years is summarized in table [S] 

Table 3. Detection efficiency e in % as a function of the 
timescale in days for events with uq < 1.5 and normal- 
ized to Tobs = 5 yr. 
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4. Blending effect 

The above efficiencies have been derived without taking 
account of stellar blending. 

We have performed a study to quantify precisely the 
influence of blending on our results. We simulate an ex- 
periment with synthetic images, regularly spaced in time, 
in which microlensing events are embedded. The images 
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Fig. 2. Red and blue light curves of the microlensing can- 
didates detected towards the SMC, in 25-day bins. 



are built from the hst colour-magnitude diagram (CMD) 
and luminosity function (LF) of the Magellanic Clouds, 
which extend well beyond the EROS detection limit (down 
to V=23.5 in the hst data). We have checked that the 
CMD and LF that we construct from these simulated im- 
ages are in good agreement with those obtained from our 
data. All stars are lensed for a timescale of 5 images. 
The impact parameter uq is drawn uniformly between 
and 1.5, and to is set so that no two stars less than 20 pixels 
apart are magnified simultaneously (typical EROS images 
have a seeing of 2.1 arcsec corresponding to a gaussian 
with rms of 1.5 pixel). 

The standard photometry chain of EROS is applied 
to these images, leading to the detection of 5599 objects 
(in this terminology an object usually encompasses several 
simulated stars) with their red and blue light curves. A to- 
tal of 6304 generated microlensing events were found with 
Uq < 1 and an average timescale {tE) — 3.16 images. The 
effect of blending is readily seen from the fact that more 
events than objects are detected. Of these events, 60% are 
due to the brightest star in the two pixels around the ob- 
ject, with recovered (^e) = 4.00 images and (uq) = 0.52, 
close to the generated values. The remaining events are 
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due to a fainter, blended component of the object, with 
recovered (^e) — 1-89 images, well underestimated, and 
{uq) = 0.64, clearly overestimated, as expected from the 
impact of a significant blend. The optical depth being pro- 
portional to the number of events and to their mean tg, 
an estimate of the ratio R of the recovered optical depth 
to the generated one is: 



R 



6304 3.16 
5599/1. 5"5ir 



= 1.07 



with an estimated error of about 10%. 

The light curves from the above blending simulation 
are then used as templates in our Monte Carlo chain de- 
signed to compute the efficiency of our analysis (section 3) 
as follows: for each object detected in the actual data, an 
object is randomly chosen from the blending simulation in 
the same region of the CMD. The light curve of the latter 
is shifted and stretched to account for the randomly se- 
lected duration and time of maximum amplification, but 
the amplitudes in each colour are of course kept as found 
in the blending simulation to reflect the impact of the un- 
derlying stellar companions of the microlensed star. This 
template is used (instead of the standard shape derived 
by Paczyhski in the absence of blending and used for the 
efficiency computation in the previous section) to add a 
microlensing event to the real (flat in most cases) light 
curve. To ensure similar photometric dispersion in the 
Monte Carlo and in the data, the rms flux deviation of 
each simulated data point from the blended microlensing 
light curve is taken to be the same as the deviation of the 
real light curve from constant brightness. The modified 
light curve goes through the same analysis cuts as in the 
standard computation of the efficiency. 

The efficiency with the blending simulation described 
above is within 10% of the efficiency computed from the 
standard pro cedure, in agreement with es timates in pre- 
vious works 1 Palanaue-Delabrouillfl Il99^ . The loss due 
to the under-estimation of the magnification of a blended 
event (the observed magnification on a blended object is 
lower than the physical magnification of the underlying 
star) is almost exactly compensated by the fact that an 
object includes more than one star subject to lensing. In 
addition, as expected, the cuts of the analysis are suffi- 
ciently loose not to reject blended microlensing events. In 
particular, there are no cuts requiring a strict achromatic- 
ity of the event. For the EROS experiment, blending has 
a small impact on the efficiency and can therefore be ne- 
glected. 

5. Limits 

In order to set limits on the contribution of dark objects 
to the halo, we use the so-called "standard" halo model 
described in SM98 as model 1, and take into account the 
efficiency of the analysis given in section |21 

For a given experiment, assuming that the halo is made 
of compact objects having a single mass M, we define 
PAiitE) as the probability distribution of expected event 



durations, taking efficiencies into account, and /-/Vm as 
the expected number of events, where / is the halo mass 
fraction with respect to a full standard halo. We construct 
a frequentist confidence level 1 — p„ considering only the 
number n of observed events as a Poisson process: 



-fNM 



i=0 

We also make an independant Kolmogorov test com- 
paring the observed event durations to the expected 
PjV,/(tE) distribution, and obtain a Kolmogorov probability 
Pt- 

A well known prescription for the confidence level of 
the combined test is 



P = PnPti^ - InpnPt) 

Our experiment has however an unknown contribution 
of background — variable stars or any other unknown phe- 
nomenon — to the observed candidates. We therefore test 
all combinations, attributing candidates either to signal 
or background, and retain, for every mass fraction /, the 
configuration that maximizes p, i.e. the most conservative 
one. This ensures that we have a given frequentist cover- 
age for all possible hypotheses. 

We then set a frequentist 95% CL limit on /, taking 
into account both the number of candidates and their du- 
ration, by finding the value of / that yields Pmax = 0.05. 

Although this prescription could be used to combine 
several experiments, it gives the same weight to all and 
does not yield useful results if the sensitivities of some 
of the experiments differ significantly. Therefore, to com- 
bine various e xperiments, we use instead a Liptak-Stouffer 
l|Lir)takLll958j) prescription where the flat p- values are con- 
verted to unit- variance, zero-mean gaussian distributed 
variables. A weighted sum of these variables is also gaus- 
sian distributed, and can be converted back to a p-value. 
We take as weights the expected number of events of each 
experiment for each lens mass^. 

The results obtained by the various phases of EROS 
(which are independent experiments) are summarized in 
tablc[21for microlensing ca ndidates towards the SMC — no 
SMC candidate in erosI (jRenault et all Il998j) . 4 candi- 
dates in EROS2 (sm98 and present work) — and in tabled] 
for microlensing candi dates towards the LMC — ii o plan- 
etary mass candidate llR,ena,u1t et alU mQ^' "l99^ LMC- 
1 from EROSl photograph ic plates An^irLetal], [1996 ) 
and the others from EROS2 llLasserre et al.ll2000HLassern , 
bOOOHMilsztain et ailioOljl . 

Figure |3| shows the 95 % exclusion limit derived from 
this method on /, the halo mass fraction, at any given 
mass M — i.e. assuming all deflectors in the halo have 
mass M — for the erosI CCD LMC and erosI CCD SMC, 
EROSl photographic plates, eros2 3-year LMC and eros2 



^ The final limit has a very weak dependance on the exact 
weighting used, provided unsensitive experiments get a very 
small weight. 
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Table 4. Summary of the microlensing candidates de- 
tected by EROS towards the lmc. The parameters are the 
same as in table |21 





Mo 




LMC-1 


0.44 


23 


LMC-3 


0.21 


44 


LMC-5 


0.59 


24 


LMC-6 


0.41 


35 


LMC-7 


0.30 


30 



5-year SMC experiments, and for the combination of all. 
We also show in the figure the limits that would be ob- 
tained with a single SMC event (SMC-1), then with no 
event at all in any of the EROS experiments, which in- 
dicates the overall sensibility of the EROS project, consid- 
ering presently analysed data. 




Mass of the deflector (Mg) 

Fig. 3. Exclusion diagram at 95 % C.L. for the standard 
halo model (4 x 10^^ Mq inside 50 kpc). The dashed lines 
are the limits towards the LMC by EROS 1 and EROS 2, 
the thin plain lines are limits towards the SMC and the 
thick line is the combined limit from the five EROS sub- 
experiments. The dotted lines are the limits that would be 
obtained considering no observed events: they indicate the 
overall sensibility of the EROS dataset. We also indicate 
the limit that would be obtained if the three very long 
duration events on the SMC were considered as background 
(label "SMC-1 only") : the limit is lowered in a negligible 
way in the 2-3 Mq mass range. 

The "dent" in the erosI plate hmit and in the eros2 
LMC limit at a mass near O.5M0 is the impact of the ~30 
day candidates observed towards the LMC. For any mass 
between 2 x 10"'^ and 10"^ Mq, we exclude at 95 % C.L. 
that more than 20% of the mass of a standard halo be 
made of compact objects. It can be seen that the com- 
bined limit is above the best limit for some values of the 
mass. This occurs quite naturally since the observed SMC 
and LMC candidates have quite different characteristics. 
At high mass, for instance, our method will consider SMC 
candidates as signal and LMC candidates as background, 
thus weakening the limit obtained by the LMC alone. It 



illustrates the (marginally significant) incompatibility be- 
tween the candidates observed by EROS towards the SMC 
and the LMC. 



6. Comparison of LMC and SMC results 

Another way to gauge the present results on microlens- 
ing towards the SMC is to compare them directly to 
the result of the MAC HO collaboration towards the 
i l200ft . 



LMC l|Alcock et al.L 1200(1) . Such a comparison is made eas- 
ier by the fact that the expected distributions of halo 
microlenses towards the SMC and LMC are very similar: 
their averages and widths differ by only a few percent, 
and these differences are much smaller than the widths 
themselves. The ratio of the optical depths towards the 
SMC and LMC is more uncertain; for a spherical isother- 
mal halo, it is expected to be about tsmc /tlmc — 1-4, 
but this value can be lower for flattened halos. Actually, it 
was proposed to use observations tow ards the SMC to eval- 
uate the galactic dark halo flattening ( Sackett and Gouldl . 

In their analysis A, the MACHO group has listed 13 
microlens candidates with an average duration < Ie > 
— 36 days and a width of the tE distribution compatible 
with all microlenses having the same mass. From this, they 
compute an optical depth t = 1.2 x 10^'', of which about 
1.0 X 10~^ is attributed to halo lenses. This translates into 
an expected number of events for the present EROS SMC 
data set of 



<tE > TT 



— < e > T SMC 



where iV* = 5.2 10^ is the number of SMC stars moni- 
tored by EROS, Tobs — 5 yrs is the total duration of the 
observations, and < e >= 0.146 is the average EROS 
detection efficiency for microlensing events of duration 
< tE >— 0.1 yr. (This average efficiency is obtained 
by interpolating the values in Table 3, and is normalized 
to uq < 1.) Substituting these numerical values leads to 
N,,ts = 2.4 X (rsMc/lO-^). 

There are, however, no EROS SMC microlensing can- 
didates that are compatible with the tE distribution ob- 
served by the MACHO group towards the LMC: candidate 
SMC-1, the shortest of all, is very likely due to a lens in 
the SMC (see parallax analysis and the discussion that 
follows in SM98); the other 3 candidates all have dura- 
tions longer than 240 days, clearly incompatible with an 
average of 36 days. From this absence of microlensing 
candidates with the required duration, EROS can exclude 
TSMC > 10^^ for events similar to those in MACHO sample 
A, at better than 90% C.L. The value corresponding to the 
spherical isothermal halo model, tsmc = 1-4 x temc — 
1.4 X 10""^ is excluded at better than 96% C.L. 

7. Discussion and conclusions 

The analysis of five years of microlensing data towards the 
SMC yielded four candidates, and allowed to put stringent 
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limits on the amount of galactic dark matter made of com- 
pact objects. Objects with a mass between 2 x 10"'' 
and 1 Mq cannot account for more than 25% of the mass 
of a standard spherical, isothermal and isotropic Galactic 
halo of 4 X 10^1 Mq out to 50 kpc. 

The new method described in this paper to combine 
the results and use all the information available on the 
events has one drawback compared to methods used in 
previous works: it could be more sensitive to the halo 
model, adding an extra dependence on the velocity dis- 
tribution. Previously, only the influence of the halo model 
on the number of expected events has been studied, as 
in SM98. A change in the velocity distribution, however, 
would mainly shift the "dent" in the EROS2 LMC limit. In 
addition, since the final exclusion limit is quite flat over a 
large range of masses, the flnal result should not be influ- 
enced by reasonable changes in the velocity distribution 
assumptions. 

This limit excludes a sizeable fraction of the allowed 
domain shown in Fig. 12 of Alcock et al. (2000). 

All the candidates that have been detected so far 
towards the SMC have long durations, and seem more 
compatible with unidentifled variable stars or self-lensing 
within the cloud than with halo objects. These would need 
to have supersolar masses to account for such durations. 

The statistics are still very low for a quantitative com- 
parison of the SMC and LMC samples of microlcnsing candi- 
dates. They have, however, quite different characteristics 
which is poorly compatible with a unique population of 
lenses. 
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